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Abstract 

We calculate the amplitude of the rare flavour-changing neutral-current de¬ 
cay B —>■ 7r£~^£~ at large recoil of the pion. The nonlocal contributions in 
which the weak effective operators are combined with the electromagnetic 
lepton-pair emission are systematically taken into account. These ampli¬ 
tudes are calculated at off-shell values of the lepton-pair mass squared, 
A < 0, employing the operator-product expansion, QCD factorization and 
light-cone sum rules. The results are htted to hadronic dispersion relations 
in including the intermediate vector meson contributions. The disper¬ 
sion relations are then used in the physical region > 0. Our main result 
is the pro cess-dependent addition AC^^\q^) to the Wilson coefficient (Tg 
obtained at 4m^ < ^ Together with the B ^ tt form factors 

from light-cone sum rules, this quantity is used to predict the differential 
rate, direct CP-asymmetry and isospin asymmetry in i? —)■ n£~^£~. We also 
estimate the total rate of the rare decay B —)■ vrz/P. 


^on leave from the Department of Theoretical Physics, Yaroslavl State University, Russia 
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1 Introduction 


The first measurement of the —)■ 7r+/i+/x“ decay by the LHCb Collaboration 

[1] paved the way for more detailed measurements of 6 —)■ di^i~ decays. These re¬ 
sults will complement the available data on 6 —)■ decays, providing new impor¬ 

tant insight in the dynamics of flavour-changing neutral-current (FCNC) transitions 
in Standard Model (SM) and beyond. 

One important feature of exclusive b —)■ d£~^£~ decays is a non-vanishing direct 
CP-asymmetry. In Standard Model (SM) this effect is caused by the interference 
between the dominant short-distance contributions of semileptonic and magnetic 
dipole operators and the contributions of other effective operators accompanied by 
the electromagnetic lepton-pair emission. The amplitudes of the latter contributions 
are process-dependent, and are defined as hadronic matrix elements of nonlocal oper¬ 
ator products. Importantly, the parts of the exclusive b —)• d£^£~ decay amplitudes 
proportional to = VubV*^ and Ac = VcbV*^ are of the same order of Cabibbo 
suppression and, in addition to a relative CKM phase, have different strong phases 
originating from the nonlocal amplitudes. The main goal of this work is to calculate 
the hadronic matrix elements of nonlocal contributions to 5 —)■ 7r£~^£~ at large recoil 
of the pion, that is, at small and intermediate lepton-pair mass, -C m% . 

An advanced theoretical description of the exclusive semileptonic FCNC decays 
was developed on the basis of QCD factorization (QCDF) approach [2], applied first 
to the decays B —)■ K^*H^£~ in Ref. [3] and to R —)■ p£'^£~ in Ref. |1]; see also further 
applications to R —)■ K£^£~ [5], and to R —)■ 7i£'^£~ [6]. An approach combining 
QCD light-cone sum rules with QCDF at > 0 for R —)■ K^*H^£~ was used in [7]. 

In QCDF, the nonlocal effects in these decays are described in terms of hard- 
scattering quark-gluon amplitudes with virtual photon emission, convoluted with 
light-cone distribution amplitudes (DAs) of the initial R meson and final light meson. 
Soft gluons, responsible for the onset of long-distance effects in the channel of the 
electromagnetic current, including vector resonance formation and nonfactorizable 
interactions with initial and final meson remain beyond the reach of QCDF. Hence, 
these contributions have to be kept small, protected by their power suppression. In 
particular, one avoids the intervals of lepton-pair mass squared g^ in the vicinity 
of vector meson masses, g^ ~ my {V = p, cu,..., J/"^,..), where nonlocal effects are 
largely influenced by long-distance quark-gluon dynamics. This constraint defines 
the region of applicability of QCDF, that is, roughly from g ^,^ = 2 GeV^ up to g^^^^ = 
6 GeV^. In this region, quark-hadron duality approximation is tacitly assumed 
for the contributions of radially excited and continuum hadronic states with the 
quantum numbers of light vector mesons. 

Note that in the R —)■ '7i£^£~ decay, as compared to R —)■ K£~^£~, the role 
of nonlocal effects related to p and u resonances in the g^ channel grows due to 
the current-current operators with large Wilson coefficients in the ~ \u part. For 
the same reason, the weak annihilation combined with virtual photon emission, 
being suppressed in R —)■ K£~^£~ decays, becomes one of the dominant nonlocal 
effects in R —)■ tt£'^£~ . In the QGDF approach one describes the weak annihilation 
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contribution [3] in terms of a virtual photon emission off the spectator antiquark 
in the B meson, followed by the subsequent annihilation to a hnal pion state. The 
accuracy of this leading-power diagram approximation, presumably quite sufficient 
for B —)■ Ki^i~ decay, becomes crucial for B —)■ 

In this paper we calculate the nonlocal effects in i? —)■ using the method 

formulated in Ref. [H] and applied in Ref. [H] to R —One avoids applying 
QCDF directly in the physical region > 0 and calculates the amplitudes of nonlo¬ 
cal contributions at deep spacelike < 0, |g^| where the operator-product 

expansion (OPE) and QCDF can safely be used. The OPE contributions include 
the leading-order (LO) loops and weak annihilation, the NLO perturbative correc¬ 
tions to the loops and hard spectator scattering. Furthermore, we include important 
nonfactorizable soft-gluon effects via dedicated LCSR calculations of hadronic ma¬ 
trix elements. The amplitudes of nonlocal effects are then represented in a form of 
hadronic dispersion relations in the variable where vector mesons are included ex¬ 
plicitly. The residues of vector-meson poles related to nonleptonic R —)■ Evr decays 
are hxed, using experimental data and/or QCDF estimates. The nonresonant part 
of the hadronic dispersion integral is parametrized combining quark-hadron dual¬ 
ity with a polynomial ansatz. Finally, the unknown parameters in the dispersion 
relation, most importantly, the strong phases of resonance and nonresonant contri¬ 
butions are htted to the QCD calculation at g^ < 0. The advantage of describing 
nonlocal contributions to the FCNC decay amplitude in terms of hadronic disper¬ 
sion relation is that the latter is valid in the whole large-recoil region specihed as 
4mj < g2 < 

In R —>■ 7Ti~^£~ decays the combinations of CKM factors Xu and Ac are compa¬ 
rable in size. Correspondingly, we have to calculate separately two hadronic matrix 
elements of nonlocal effects multiplying and Ac. A similar CKM separation has 
to be done in the amplitudes of nonleptonic R —)■ Dvr decays, used to £x the residues 
of vector-meson poles in the hadronic dispersion relations. To obtain the separate 
parts of nonleptonic R — )■ Rtt amplitudes for V = p,u we employ the QCDF results 
[To] and control the resulting amplitudes with the data on branching ratios and 
CP-asymmetries of these nonleptonic decays. 

The plan of this paper is as follows. In Sec. we present the structure of the 
R —)■ 'k£'^£~ decay amplitude and dehne the hadronic matrix element of nonlocal 
contributions. Sec. contains a detailed calculation of these amplitudes at g^ < 0. 
In Sec. I^we perform the relevant numerical analysis. In Sec. [^the necessary inputs 
for the nonleptonic R —)■ Dvr decay amplitudes are presented. Sec. |^is devoted to 
the analysis of the hadronic dispersion relations. Matching the latter to the result of 
QCD calculation, we then obtain ACQ^^\q‘^ >0). In Sec. [^our predictions for the 
observables in the R —)■ decay are presented, including the decay rate, direct 

CR-asymmetry and isospin-asymmetry. In Sec. [^we estimate the rate for R —)■ ttz/P 
decay and Sec. contains the concluding discussion. The two appendices contain: 
(A) the operators and Wilson coefficients of the effective Hamiltonian of 6 —)■ d£^i~ 
transitions and (B) the QCDF expressions used for the amplitudes of R —)■ p(ci;)7r 
nonleptonic decays. 


3 





Figure 1: FCNC contributions to B ^ due to the effective operators Og^io 

(left) and Oj-y (right) denoted as black squares. 


2 The B decay amplitude 


The effective weak Hamiltonian of the b —)■ d£^i transitions {£ = e, /r, r) has 
the following form mm in the SM : 


Tjb^d 

-^eff 


AGf 
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+ ]+h.c. 


( 1 ) 


2=1 


2 = 1 


2=3 


where Xp = Vp^Vf^, (p = u, c, t) are the products of CKM matrix elements. In 
contrast to the b —)■ s£^£~ transitions, all three terms in the unitary relation have 
the same order of Cabibbo suppression, A^ ~ Ac ~ At ~ A^, A being the Wolfenstein 
parameter. Hereafter, we assume CKM unitarity and replace At = — (A^ + Ac). The 
local dimension-6 operators Oi in Q together with the numerical values of their 
Wilson coefficients Ct at relevant scales are presented in the Appendix A. 

The amplitude of the B —)■ decay reads: 


A{B^7i£+£-) = -{7i{p)£+£-\H^^ff'^\B{p + q)) 

GF Clem 


TT 


-At 


(h‘‘e)pAc,!M) + 


2mb 


CffLiQ^) 


tub + ruj, 


+ 


(<V75<) ) 


.( 2 ) 


where p^ and are the four-momenta of the vr-meson and lepton pair, respectively. 

In the dominant contributions of the operators Og^o and Ojy are separated 
(Fig. 1^ and their hadronic matrix elements are expressed in terms of the vector 
and tensor H vr form factors, and /J.n-(5'^), respectively, dehned in the 


4 






















standard way: 


{■K{p)\dYh\B{p + q)) 
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ml 
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( 3 ) 


{Ti{p)\da^'' q^h\B{p + q)) 


mB + 


2q^p^^ + 



(m| - ml) 



(4) 


For definiteness, hereafter we consider the B~ —)■ mode, nnless stated oth¬ 

erwise. We assnme isospin symmetry for the b ^ d and b ^ u transition form 
factors. The B~ —)• tt~ form factor in Eq. ([^ is eqnal to the one in the 
B^ —> semileptonic decay and the form factors in the B^ —)• decay 

amplitnde have an extra factor 1/\/2. For the CP-conjugated modes B^ —)■ vr+f'+f'" 
and P° —)■ , respectively, one has to use the hermitian conjugated effective 

operators with complex conjugated CKM factors A*. 

The current-current, quark-penguin and chromomagnetic operators in the ef¬ 
fective Hamiltonian Q contribute to the decay amplitude (|^, with the lepton 
pair produced via virtual photon. After factorizing out the lepton pair, the ex¬ 
pression for these nonlocal effects is arranged in ([^ in a form of correlation func¬ 
tions of the time-ordered product of quark operators with the quark e.m. current, 
i/i™ = J2q=udscbQQ^'y^J'^y saudwiched between B and vr states: 




i / d^xP‘^^(7r(p)|T 



CiCj’(O) + C2Ci(0) 


+ CfcCfc(O) 

k=3—6,8g 


\B{p + q)) = [{p-q)qf. 




u,c), (5) 


where the index p = u,c hereafter distinguishes the hadronic matrix elements in 
Eq.(j^, multiplying, respectively, the CKM factors Xu, K- Substituting (|^ in (|^ 
and taking into account the conservation of the leptonic current, we write down the 
decay amplitude in a more compact form: 




(ere) pJc, + Acf’>(«") 


2mb 


ms + m, 


+ {h^ld)pgCiQ 


( 6 ) 


where the invariant amplitudes introduced in Eq. ([^ form a process-dependent and 
g^-dependent addition to the Wilson coefficient Cg: 


ACf = -I67r2 


{XuU^^Kq^) + XcU^^Kq^)) 


(7) 
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Figure 2: Leading-order diagrams of nonlocal effects in B ^ due to the four- 

guark effective operators Off O^-q: the guark-loop (a) and weak annihilation 
(b). The black sguare denotes the operator and crossed circles indicate possible points 
of the virtual photon emission. 


In Eq.([^ we also introduce the ratio of tensor and vector form factors: 

rLif) - ■ ( 8 ) 

fBnir) 

In the case of 6 —)■ transitions, the factor is usually neglected so that Ac = 

—Xt, and one recovers the corresponding expression for ACg^^\q‘^) in B ^ Ki^£~ 
used in Ref. j9|. 


3 Nonlocal effects at spacelike 

In this section we present separate contributions to the nonlocal amplitudes 
and dehned in ([^ and calculated at < 0, in the same approxi¬ 
mation that was adopted in Ref. [9] for B —)• . 


3.1 Factorizable loops 

At leading order (LO) in the quark-gluon coupling, the contributions of the 
four-quark operators to R —)■ have two possible quark topologies. One of 

them corresponds to the factorizable quark-loop diagrams with different flavours 
(see Fig. [^a)). Their expressions are obtained from, e.g., the ones presented in 
Ref. [9], separating the ~ A„ and ~ Ac parts: 


{u 
fact,LO 




-i/ic) 

^fact,LO 




9o{(f)fBM^) + > ( 9 ) 

^ Y + giq\ rnDf+^iq"^) + nll~%{q ^), (10) 
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where the common term stemming from the quark-penguin operators Os-e is 


a /( 3 - 6 ) (^2 
^fact.LO 


= 


247r2 


+ ^^'4 + C5 - 1 - ^ '^D + ■, '^Vi) 


+ 2 ^Cs -f -(^4 -I- C 5 -1- ml) 

+ ^C's + -^ + C 5 -f 9o{q^) 

C's + -^ + C's + ■ 


( 11 ) 


For the loop function we use the expression valid at < 0: 
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X In 
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+ 3-‘V^ 


- a/^ 
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1 _ 




2ml 

—^ + 1 
qz 


( 12 ) 


where is the quark mass if g = 6, c, s and /i is the renormalization scale. For 
u- and d-quark loops the quark masses are neglected; in this case the loop function 
takes the form: 


9 o{q ) = lim 9(9 , ml) = - - In 

mj —^0 O 


-q 




(13) 


In Eqs. and (10) the “full” I? —)• vr form factor is the same as in the contributions 


of Og^io operators to the decay amplitude (|^. For this form factor we will use LCSR 
results that are valid also at < 0. 

Note that in the LO approximation, when gluon exchanges between the loop 
and the rest of the diagram in Fig |^a) are neglected, the nonlocal amplitudes 
^iacfhoiq'^) I*® calculated within LCSR approach. One has to dehne the 

vacuum-to-pion 3-point correlation function of the R-meson interpolating current, 
the four-quark operator and the electromagnetic current. After the quark loop is 
factorized out at large spacelike g^, the remaining correlation function coincides with 
the one used to calculate the R —)■ tt form factor from LCSR. The resulting sum 
rule is then reduced to the loop factor multiplied by the LCSR expression for the 
R —)■ TT form factor, reproducing Eqs. §, Q. 


3.2 Weak annihilation 

The second possible topology at LO is the weak annihilation (WA) with the 
diagrams shown in Fig. |^b). In QCDF, neglecting the inverse heavy 6-quark mass 
corrections, the leading diagram is the one where the virtual photon is emitted off 
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the spectator quark q = u^d m. the B meson, with the resulting expression HE]: 




1 fBUrrib 
8Nc 



J du cu), 

0 


(p 


u,c), 


(14) 


where and fs are the vr- and i?-meson decay constants, respectively, and the 
hard-scattering amplitude 


t1°)’^(u,cu) = Q,C'P 


ArriB tubu: 

WA 9 

rUb tubOO — 


(P = M, C) 


(15) 


is convoluted with the 5-meson DA dehned as in Refs. [131 El and (Ptt{u) is 

the twist-2 pion DA. The factor 


C^A = ^u{SUC2 + 3Ci) + S^Ci + 3 C 2 )) + Cs + 3^4 (16) 


is the combination of Wilson coefficients depending on the ffavour-content of the 
B meson. To obtain the amplitudes 'H^p^{q‘^), one takes into account that the LO 
kernel (15) is independent of the variable u, hence the integral over (pTr{u) is reduced 
to its unit normalization. Adopting the exponential ansatz [T3| for the 5-meson 
DAs: 


4>U<^) = , (17) 

Ab 

where As is the inverse moment, we obtain the following expression for the amplitude 
valid at < 0: 


Wwa(«") 


2NcmB)^B 


mB^Bj 


r'P 

'^WA 5 


(18) 


where Ei(x) = — f dte ^/t. 

— X 

In contrast to 5 —)■ Ki^i~ transitions, the WA mechanism due to the enhanced 
current-current operators 0 ^ 2 ^ provides one of the dominant contribution to the 

^wa(^^) amplitude in 5 —)■ Moreover, the resulting difference between the 

WA amplitudes in B~ —)■ and 5° —)■ contributes to the isospin 

asymmetry in 5 —)■ . 

Since the role of WA effects becomes important, it is desirable to improve the 
accuracy beyond the leading diagram contribution considered here. We checked that 
adding all subleading diagrams in Fig. [^b) to the virtual photon emission from the 
spectator quark does not produce a visible effect for the Oi ,2 contributions. There 
still remain power suppressed corrections generated by the higher twists in the pion 
DAs, and the contributions of the operators q yet unaccounted in QCDF. In the 
future also the perturbative nonfactorizable corrections to the diagrams in Fig. |^b) 
have to be calculated. 







X 



Figure 3: Factorizable NLO quark-loop contributions to B ^ , due to the 

four-quark effective operators Of 2 (upper panel) and the chromomagnetic operator 
Ogg (lower panel). The notation is the same as in Fig. 2. 


In principle, it is also possible to calculate the WA contribution employing the 
LCSR approach with the R-meson DAs. The correlation function will be described 
by the diagram similar to Fig. Sb). but with the on-shell pion replaced by the 
interpolating quark current with the virtuality After employing the hadronic 
dispersion relation and quark-hadron duality in the pion channel, in the factorizable 
approximation, the two-point part of this correlation function will yield the QCD 
sum rule for the pion decay constant squared. The result for then yield 

the expression (18). 


3.3 Factorizable NLO contributions 

The NLO corrections to the quark loops generated by the current-current op¬ 
erators Off are given by the two-loop diagrams shown in Fig. [^( a),(b),(c). The 
contributions of quark-penguin operators are neglected, being suppressed by small 
Wilson coefficients and as simultaneously. In the same order of the perturbative 
expansion, the chromomagnetic operator Ogg is described with the diagrams shown 
in Fig. [^d),(e). These factorizable NLO contributions were taken into account in 
QCDF [31 0], employing the quark-level two-loop diagrams calculated in Ref. [13] 
(see also Refs. [lain]). The NLO contribution of Of 2 to can be literally taken 
from Ref. [H], replacing by The corresponding contribution to 

has the same structure, so that we can present both contribution by one compact 
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expression: 


njip) 

^fact,NLO 


3271^ tub 




eff t7'(7) / 2\ 


ms 

2mb 


CiFf + 2C2 



(19) 


where p = u,c. The dehnitions and nomenclature of the indices of the functions 
-p 2 ,p^^ are the same as in Refs. [T^ [T6]. The only difference is that 

and F^'^f'^ are expressed as a double expansion in s = /ml and ml = ml/ml, 
whereas f[^^^ and are expanded only in powers s = /m^ since we work in 

the limit m^ = 0. It has been shown in Refs. [T^ US] that keeping the terms up to 
the third power of s and ml provides a sufficient numerical accuracy in the region 
0.05 < s < 0.25. Here we use this expansion for q^ < 0, restricting ourselves to 
1.0 < \q^\ < 4.0 GeV^, i.e. staying within the same region. For we use the 
expression derived in Ref. |3]. 

We remind that at NLO, the nonlocal contributions acquire the imaginary part 
also at q^ < 0, that is, not related to the singularities in the variable q^. The origin 
of this imaginary part and its relation to the final-state strong interaction is the 
same as for B —)■ Ki^i~ and is explained in detail in Ref. j9|. 

Note that analytic expressions for the two-loop virtual corrections to the matrix 
elements of the 0“ and O 2 operators are available from Ref. Id- These expressions 
are valid at > 0 and agree with the expansion in q^/ml obtained in Ref. jl6j . 
However, we cannot use the results of Ref. na straightforwardly in our calculation 
at < 0, without separating the imaginary contributions inherent to the negative 
g2-region from the contributions appearing due to the cuts of quark-gluon diagrams 
at q^ > 0. Hence, we prefer to use the expanded form of these corrections [T6| in 
which the phases stemming from the positivity of g^, e.g., the terms proportional to 
ZTT and originating from the log g^ terms can be easily recognized and separated. As 
we work at sufficiently small values of | g^ |, the accuracy of the expansion in Ref. [16] 
is sufficient. 

Contrary to the LO contributions considered in the previous subsections, the fac- 
torizable NLO ones are not simply accessible within the LCSR approach. Indeed, in 
order to reach the same 0 {as) accuracy, the calculation of the underlying correlation 
function has to include two-loop diagrams with several scales, a task exceeding the 
currently reached level of complexity in the multiloop calculations. 


3.4 Nonfactorizable soft-gluon contributions 

We also take into account the nonfactorizable contributions to the amplitudes 
'^(p)(g2) emerging due to a soft-gluon emission from the quark loops, as shown 
in Fig. These hadronic matrix elements cannot be reduced to a R —?• vr form 
factor. It is also not possible to attribute the soft gluon to one of the hadrons in 
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Figure 4: Nonfactorizable soft-gluon contributions to B ^ due to (a) the 

four-quark effective operators Of 2 and O 3-6 and (b) the chromomagnetic Ogg op¬ 
erator. The soft gluon is represented by the gluon line with a cross. The rest of 
notation is the same as in the previous figures. 


the B ^ 71 transition. The soft-gluon contributions are nevertheless well dehned at 
< 0 and |g^| 3> shown in Ref. [8], their suppression with respect to the 

factorizable loops is controlled by the powers of l/(4m^ — q^) and l/|g^|, stemming, 
respectively, from the c-quark and massless loops with soft gluon. The corresponding 
hadronic matrix elements were hrst calculated for the c-quark loops in Ref. [8] for 
B —)■ The calculation was done in two steps: (1) applying the light- 

cone OPE at deep spacelike for the quark loop with soft-gluon emission, and (2) 
calculating the B —)■ K^*'> hadronic matrix element of the emerging quark-antiquark- 
gluon operator using LCSRs with the R-meson three-particle DAs. Completing this 
result to include the loops with all possible quark flavours is straightforward and was 
already done for B —)■ Ki^i~ in Ref. [9]; to obtain the corresponding contribution 
to 'H^^\q^) in R —)■ , we only need to replace the kaon by the pion. The 

soft-gluon nonfactorizable contribution of the operator Of contributing to 'H^'^\q‘^), 
is also easily obtained. The result for this contribution is cast in a compact form: 

= fii^C,+Ci-Cs)A{mlq^) + fi25,„C,+Ci-Cs)A{0,q^) 

- fic3 + C4-Ce)(^A(m^„q‘^) + A(ml,qfi. ( 20 ) 


A cumbersome expression for the nonfactorizable hadronic matrix element A(mg, g^) 
obtained from LCSR can be found in Ref. [S],(see Eq. (4.8) therein), where the 
dependence on the quark mass squared is explicitly shown and is indicated in the 
above expression. To adjust this expression to the R —)■ transition, one has to 

replace the decay constant, meson mass and threshold parameter in this equation: 
fx —t /jT, mx —t Sq —)■ Sq, thus taking into account the flavour SU{3) violation. 

In the sum rule for g^), we use the ansatz for the three-particle R-meson DAs 


suggested in Ref. |23], with the parameter directly related to the 

inverse moment of the two-particle DA specified in Eq. (17). 
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Figure 5: Nonfactorizable spectator contributions to B ^ due to (a) the 

four-quark effective operators Of 2 , Oz-% and (b) the chromomagnetic operator O^g- 


The soft-gluon contribution of the chromomagnetic operator O^g is described by 
the diagram in Fig. |^, where instead of the loop factor, one has a pointlike emission 
of the soft gluon held. One modihes the correlation function accordingly and arrives 
at the LCSR that was already derived in Ref. [H] and presented in Eq. (4.7) therein]^ 
Making the necessary replacements for B —)■ 7 ri^i~, we obtain 


K 


{u) 

softjOsg 


(g^) = n 


(c) 

softjOsg 




^(BK) 

^softjOgg 





( 21 ) 


As in the case of R —)■ Kt^^ transitions, this contribution turns out to be very 

small. 


3.5 Nonfactorizable spectator scattering 


An important nonlocal contribution to the B —)• amplitude in NLO 

emerges due to a hard gluon emitted from the intermediate quark loop or from 
the Osg-operator vertex, and absorbed by the spectator quark in the R —)■ vr tran¬ 
sition, as shown in Fig. Following [9], we will use the QCDF result |3] for this 
contribution. The following expression is valid for both p = u and p = c parts of 
the nonlocal amplitude: 


^nonfjSpect v” / 


+ 


OisCp fBfnnVb 
327rA'c nr'^B 



j du(pT,{u)T^'*’^{u,uj)-\- 
0 



1 


/ 


du (pTr{u) uj) 



( 22 ) 


The hard-scattering kernels entering the above expression have the form: 




{u,u) 


+ Q - Ce) - it||(u,m,)(C'3 + C,- C^) - 

mb \o 3 

^■ 5 ) (^3 04 — Cg) -l- -(2SpuCi -h C 4 — (76)^11 (u, 0) I, (23) 


^We notice that in the related Eq. (4.4) a factor Csg on the r.h.s. is missing. 
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= -Q 


niBUJ 


^niBOJ — — ie 


SmB 

Srrib 


h{Tn^,UTnB uq )(5pcC*i “1“ C*4 “1“ C'e) “1“ 


-\-h{Tn ^, uTTi ^ + uq ) {C^ + C *4 + C*g) + 

+/i(0, UTn?Q + ((5puCi + (^3 + 3 C 4 + 3C*g) — 

- a - 15^6)) + ^ , 

3 J u + uq^/m^^ 


(24) 


where Qg is the electric charge of the spectator quark in the i?-meson (g = u, d) and 
the functions t\\{u,mq) and h{rn?q,q^) can be found in Ref. ^ 

The two-particle R-meson DAs 4)%{u:) are given in Eq.(17); for the twist-2 pion 
DA we employ the standard Gegenbauer expansion: 

= 6m(1 - u)(^ + a2{iji)cf^‘^\u) + al{n)cf^‘^\u)^ . (25) 


The fact that the amplitudes in (22) depend on the charge of the spectator quark 


in the B meson triggers another important contribution to the isospin asymmetry 
in R —)■ . 

Summing up all contributions considered in this section, we obtain the two non¬ 
local amplitudes in the adopted approximation: 


= W[h,Lo(9' 


fact,LOVy ) T‘H^^(g j -r ' i-fact.NLO 


+ W2.(9t + C!.,0.(9t + W”,,p=c,(9"). (P = n,c). (26) 


/(P) 


dp) 

Uac 

/(P) 


q-lVP) (n'^\ 


4 Numerical analysis 


Here we perform the numerical analysis of the nonlocal amplitudes (26) at space¬ 
like q^ < 0, more dehnitely, in the region 1 GeV^ < |g^| ^ 4 GeV^, where the OPE 
and QGDF approximation can be trusted. The input parameters and references 
to their source are listed in Table the charged and neutral R-meson and pion 
masses are taken from [TB], and the numerical values of the Wilson coefficients are 
presented in the Appendix A. As a default renormalization and factorization scale 
we assume p = 3 GeV, the same as in Ref. [9]. It will be varied in the interval 
2.5 < p < 4.5 GeV to study the fi dependence. 

For the vector R —)■ vr form factor the most recent update [19] of LGSR prediction 
is adopted, in a form fitted to the three-parameter BGL parameterization: 




1 — q^ lm\ 


B* 


z{q^,to) - z(0, to) - ^ (^z{q‘^, tg)^ - ^(0, fg)^) 


+bt 


z{q‘^, to)^ - ^(0, to)^ + ^ (^(g^ to)^ - ^(0, tg)^) 


(27) 
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Figure 6: Nonlocal amplitude for B —>■ vr at < 0 calculated at the 

central values of the input; the real (imaginary) part is in the upper (lower) panel; 
^he sum of the separate contributions specified in Eq.(26). 


with to = (u7.B + m,r)(\/’^s~\/^^ 7 r)^t where the normalization and shape parameters 
are presented in Table The decay constant of B meson is determined from two- 
point sum rules, where we use the recent analysis izni; the inverse moment of the 
5-meson DA is also represented by the interval of QCD sum rule prediction |21j . 
The intervals of Gegenbauer moments in the pion DAs used in the QCDF expressions 
are the same as in the LCSR for the 5 —vr form factor mm- 

Substituting the central input in the expressions presented in the previous sec¬ 
tion, we calculate the two amplitudes < 0) and 'H'^^\q^ < 0) for B~ —)■ 

and plot their real and imaginary parts in Figs. |^and[^ respectively, show¬ 
ing also the separate contributions. For comparison, the amplitude < 0) is 

plotted in Fig. ^for B^ —)■ whereas the amplitude < 0) for the latter 

mode is numerically similar to the one for B~ —)■ and is not shown. 
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(GeV2) 



- Im[HS] 


Figure 7: The same as in Fig'll /or the amplitude for B —)■ vr £+£ 


From the numerical analysis we can draw several conclusions: 

(a) The contributions to TL^‘^\q^) are approximately the same as the correspond¬ 
ing ones for the B —)■ K£~^£~ obtained in Ref. |9j, the differences reflect the violation 
of the flavour SU (3) symmetry; 

(b) The contributions to 'H^^\q‘^), in B -^71 £^£ presented in Fig. ^ are 
clearly dominated by the weak annihilation term, enhancing the real part of TL^{q^) 
considerably. This effect is less pronounced for —)■ 7 i~£^£~, as expected; 

(c) The nonfactorizable soft-gluon contribution due to the operators 0“2 are im¬ 
portant and the corresponding contributions of O 3-6 are not negligible. Meanwhile, 
the contribution due to the operator O^g with a soft-gluon is very small. 

The uncertainties of the functions TL^^\q^ < 0) and 77^“^< 0) are estimated 
varying the inputs within their adopted intervals indicated in Table The largest 
uncertainties originate from the variation of Jb, As and the correlated variation of 
the parameters of To stay on the conservative side, we neglect possible corre¬ 
lations between the individual input entries in Table [1} We also varied the renor- 
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Figure 8: 


The same as in 


Fig^for the amplitude 


in —)■ . 


malization/factorization scale around the default “optimal” value p = 3.0 GeV. The 
results do not signihcantly change the estimated total uncertainty and we therefore 
neglect the scale dependence in the error estimates performed below. 


5 Nonleptonic B ^ Vtt decay amplitudes 

We now turn to weak nonleptonic B ^ Vir decays with neutral vector mesons 
V = p^,oj,(j), J/'ip,'ip{2S) in the hnal state. The intervals of the absolute values 
of their amplitudes have to be estimated and used as an input in the hadronic 
dispersion relations for to be htted to the calculated < g). 

The amplitude of a B~ V7i~ decay is parametrized as: 

AiB-^VTT-) = {V{q)7r-ip)\H^J^^^^\B-{p + q)) 

= ^ "^vi^v ■ p) ) ( 28 ) 
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Input parameter 

[Ref.] 

a,{mz) = 0.1185 ±0.0006 

mc(mc) = (1-275 ± 0.025) GeV 
rhbijnb) = (4.18 ± 0.03) GeV 
m,(2GeV) = (95 ± 5) MeV 

jllBI 

fB = 2071^" MeV 

m 

Xb = (460 ± 110) MeV 


u = 130.4 MeV 

ra 

aJ(lGeV) = 0.17 ±0.08 

<(lGeV) = 0.06 ±0.10 

ra 

Sum rules in the pion channel 

M2 = 1.0 ± 0.5 GeV2 , = 0.7 GeV^ 


/+^(0) = 0.307 ±0.020 
bf = -1.31 ±0.42 
b+ = -0.904 ± 0.444 

jra 


Table 1: Intervals of the input parameters used in the calculation of < 0). 


where q {ey) is the 4-momentum (polarization vector) of the vector meson with = 
my. For the charge-conjugated mode Vti^ one has to replace 


K,c ill 

U^C 


the above relation, whereas the hadronic amplitudes remain unchanged, = 

A 


U,C 

B+V7r+' 


amplitudes as 


mg 


For the neutral —?• decay modes we denote the corresponc 

AA^Qy^Q. The effective Hamiltonian i'l 

contains the operators O 3 - 6 ) O^g given in Appendix A, and we neglect 

electroweak quark-penguin operators with 0(Q!em) suppressed Wilson coefficients. 
From Eq. (28) one obtains the expression for the CP-averaged width; 


(28) 

the 


r(B^ [r(B- ^Vn-) + r(B+ ^ >4. mp 

Zi oTT/TT-^ 

X (^^u\‘^\A'^-y^-\‘^ + lAcPlA^-y^-p -|- 2|AuAc||A^-y^-||A^_y^_| cos(5cos , (29) 


where A (a, b, c) = + A — 2ab — 2ac — 2bc. In the above, we explicitly isolate 

the relative strong phase 6 between the amplitudes A'f_y^_ and -A^-y^- and denote 
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A = arg(A„) — arg(Ac). The direct CP asymmetry takes the following form: 

V{B- Vtt-) - T{B+ l/7r+) 


Acp{B^ Vtt^) = 


= —2 sin (5 sin AI 


( 

\ Au 

L 

^cA^B-Vtt- 

1 

^cA^B-Vtt- 


\ Au 

AuA±b-V7V- 


2T{B^ Vtt^) 

+ 2 cos S cos A 


-1 


(30) 

Analogous expressions are obtained for the neutral B ^ Vtt^ modes, replacing 
B- B^ and B+ 5°. 

For the input in the two dispersion relations for the amplitudes and 

we need to separately determine the moduli of the hadronic amplitudes 


-Vtt 


_| and |A^_ 


B-Vtt- 


In principle, one can use the two observables in Eqs (28) 
and (30), but the presence of the third unknown parameter, the relative strong 
phase, hinders the determination. Therefore, the situation is more complex here 
than for the nonleptonic B ^ VK decays used in the analysis of i? —?■ Ki^£~ 
in Ref. [9] where only the contributions proportional to Ac were retained and the 
amplitudes were directly obtained from the measured B ^ VK branching 

fractions. 

On the other hand, there is a possibility to estimate separate contributions to 
the nonleptonic amplitudes applying the QCDF approach [2]. The latter is known 
to provide a reasonably good description of the charmless channels B~ —)■ and 

B~ —)■ . Here we use the results of Ref. |T0] where the QCDF description for 

B ^ VP decays was elaborated in detail. The necessary expressions for the am¬ 
plitude decomposition and the additional input parameters including the B ^ V 
form factors, decay constants and the Gegenbauer moments of the DAs of R = P,oj 
are collected in Appendix B. The resulting absolute values of the amplitudes are 
presented in Table To check the validity of these estimates, we calculated the 
observables (28) and (|30|, and compared the results with the experiment and with 


the earlier predictions of Ref. BDl (see Table [^, observing a reasonable agreement. 
In the transition from the widths to branching fractions we use the lifetimes of B 
mesons from Ref. [18]. Note that in the dispersion relations we will not isolate the 


Mode 

\^^Vtt\ 

\^%Vtt\ 

Mode 

\-^^Vtt\ 

\^'bVtt\ 

B^ p^TT^ 



R° —)■ 

9.9lJi 

0 

B^ —)■ UlTT^ 


n 

U-'J-O.l 


0 

0 

B^ —)■ 

n 

29.2+};^ 

R° ^ J/07r° 

D S’*"®'® 

2O.60J;? 

B^ -)■ 0(2S')7r"F 

3 

32.3^1? 

RO ^ 0(2^)7rO 

2 4+"^-’^ 

^•^-2.4 

22.8^11 


Table 2: Inputs for the absolute values of the B ^ Vtt amplitudes (in 

MeV). 


intermediate 0-meson pole, hence we do not specify the R —)• 07r nonleptonic ampli¬ 
tudes here. These decays originate either due to the q = s part of the quark-penguin 
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operators Oa-e with suppressed Wilson coefficients, or due to the Oi _2 operators 
combined with a transition via intermediate gluons into ss state. The latter is OZI 
suppressed (cf. the smallness of (j) decays into pions). The measured upper limit 
BR(i?“ — )■ (j)7r~) < 1.5 X 10“^ [18], being signihcantly smaller than the measured 
branching fractions of B~ —)■ p{u) 7 r~ decays (see Table convinces us that the 
intermediate (/)-meson contribution to R —)■ is small. Furthermore we do not 

separate the radial excitations p', approximating their contributions to the 

hadronic spectral density by the quark-hadron duality ansatz; hence we do not need 
to consider here the nonleptonic decays of the type B —)■ p'(1450)7r. 


Channel 

Observable 

Experiment 

QCDF [lO] 

QCDF, this work 

B^ 

BR X 10® 

8.3 ± 1.2 

11 . 9 ;^:? 

q c;+2.9 


Acp 

0.18^®:?? 

0.04 ±0.19 

0.09 ±0.17 

B^ 

BR X 10® 

6.9 ±0.5 

8.8^®:^ 

8.91?:® 


Acp 

-0.04 ±0.06 

-0.02 ±0.04 

-0.06 ±0.06 

BO pO^o 

BR X 10® 

2.0 ±0.5 

0 

'^•^-0.4 

0 2+®-'^ 


Acp 

— 


n 24+0-36 

RO ^ a;°7r° 

BR X 10® 

< 0.5 

n ni +0-04 
U-Ui_o.01 

n qi+ 0.06 
u-ur_o,oi 


Acp 

— 

— 

00 0 
0 d 
+ 1 

0 

1 


Table 3: Comparison of the experimental data and theoretical predictions for 
the observables in the nonleptonic B —)■ p(ci;)7r decays. 


For the neutral B^ ^ Vtt^ modes, the QCDF prediction [101 fails to predict the 
partial width B^ —)■ p^n^, the experimental value being significantly larger. Without 
going into more detailed discussion of this problem, guided by the hierarchy of 
amplitudes in the charged mode, A'^OpO.j^o ^ ^^OpO^o, we simply assume R^OpO^o = 0 
and extract I^^OpO^ol from the measured partial width employing Eq. (29). For the 
B^ —)■ uTi^ mode only the upper limit on the branching fraction is available |18j . 
indicating that this decay amplitude is suppressed in comparison to the other modes, 
hence we put ^ ~ 0 as it specihed in Table 

Turning to the charmonium channels B — >■ fjTi, where = J/'ijj,'ijj{2S), we do not 
expect the QCDF approach to work there due to a heavy hnal state and enhanced 
nonfactorizable, power suppressed effects (see, e.g., the discussion in Ref. [ 8 ]). On 
the other hand, one anticipates that these nonleptonic decays are dominated by 
the emission topology due to the operators OJ 2 with large Wilson coefficients and 
a small admixture of O 3 -Q (g = c). The contributions of the operators 0^2 
O 3 -Q (g 7 ^ c) are expected to be strongly suppressed. Theoretical estimates for 
the analogous contributions to R —)■ 'tpK transitions (see, e.g., [2l| and [25]) yield 
the amplitudes of gluonic transitions of light-quark loops to charmonium states at 
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the level of 10“^ of the dominant contributions of ^ operators. With an extra 
Cabibbo enhancement of the A„ terms in S —)■ with respect to i? —)■ 'ipK, still a 
considerable suppression remains. Hence we expect that <C |y4^_^^_|. In 

this situation the relative strong phase does not considerably influence the extraction 
of the large ~ Ac term , whereas the uncertainty of the small ~ A„ term is tolerable. 
Therefore, we use the current experimental data on the branching fractions and 
CP-asymmetries of the above decays [IB] and perform the £t of these data to the 


Eqs. (29) and (30), extracting the absolute values of the amplitudes and 

and allowing the relative phase 6 to change from 0 to 27r. The resulting 
intervals are presented in Table Finally, for the neutral P° —)■ modes, 

we make use of the isospin symmetry relation: ~ l/\/2^^-^^-, since for 

the dominant '0 emission mechanism of these decays there is only one independent 
isospin amplitude. This assumption is supported by the measurement [T8| yielding 
r(P° —)■ Jj'lpTT^) ~ l/2r(P+ —)■ J/ipTi^) . The resulting estimates are also presented 
in Table IH 


6 Hadronic dispersion relations 


Following Refs. 0191, the invariant amplitudes and are repre¬ 

sented in a form of hadronic dispersion relations in the variable inserting the 
total set of hadronic intermediate states between the electromagnetic current and 
the effective operators in the correlation functions (|^: 




(9= - 4) 


ds 




_ 1772 '^ 

y=p,a;,J/V>,P(25) ^ ^ 

A'b) 


(s - q^){s -q^ - k) 


Sh 


kyfvA^BV-K _ 

- (P - irnyV^v^) 

(p = M,c), (31) 


where the ground-state vector mesons (except 0) are isolated and the integral 
describes the contribution of excited and continuum contributions starting from 
Sh = 4m^, the lowest hadronic threshold To achieve a better convergence, we 
implement one subtraction at = —1.0 GeV^ in Fq. (31). In the above, the masses 
and total widths of the vector mesons V = cu, J/'ip,'ip{2S) are taken from Ref. 
Their decay constants are dehned as 


{0\f'^'^\V{q)) = kvmve^y{q)fv, (32) 

where the coefficients ky are determined by the valence-quark content of V and the 
quark charges: kp = l/\/2, k^ = l/(3-\/2) and = A:^( 2 S) = 2/3. The numerical 
values of fy are fixed from the measured [H] leptonic widths r(I/ —)■ £+£“) yielding 

^ Note however that a part of the 2-pion continuum contribution in this region is effectively 
absorbed in the p meson total width (for more details see, e.g., |1S|). 
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fp = 221 MeV, /(^ = 195 MeV, fj/^ = 416 MeV and /^( 2 S) = 297 MeV. The absolute 
values of the amplitudes and A^^y^ obtained from the analysis of nonleptonic 

decays in the previous section are taken from Table 

At < 0, more specihcally, in the region —4.0 GeV^ < < —1.0 GeV^, we 


substitute in the l.h.s of the relations (31) the result for specihed in Eq. (26), 

calculating simultaneously the subtraction terms at = —1.0 GeV^. The task is 
then to £t the free parameters on the r.h.s. of the hadronic dispersion relations. 
Importantly, each G-pole residue in Eq. (31) for p = u or p = c has a relative 


phase with respect to the other vector-meson contributions and to the integral over 
These phases should match the imaginary part of the calculated l.h.s. of 
the dispersion relation. As explained in Ref. [U], the phases emerge due to the 
intermediate on-shell hadronic states in the variable (p + g)^ = and are not 
related to the analytical continuation in the variable g^. 

Note that the relative strong phase between the amplitudes A'py^ 

contributing to the B ^ Vn nonleptonic amplitude, although of the same origin, is 


a different quantity, because in each of dispersion relations (31) only one of these 


amplitudes enter. On the other hand, calculating the phases of nonleptonic ampli¬ 
tudes within a theoretical framework, such as QGDF, it is possible to estimate the 
relative phase between, say, and 

In what follows, we attribute a phase to each V -pole term: 


To reduce the number of free parameters, we hx the phase differences: 






'^B~p^7r~ ^B~U7r 


_ = 0.033, 


_ X'' 

^ B~ p^7T~ ^B~Lijn~ 


dc) 


= -3.65, 


(33) 


(34) 


calculating it from QGDF, as explained in the previous section. Note that for the 
neutral mode the contribution of the w-meson is neglected and the corresponding 
difference is irrelevant. The three remaining phases and ^^( 2 S) each 

p = u,c are included into the set of £t parameters. This set will be completed below 


by the fit parameters of the integrals over p^\s). Furthermore, we adopt the Breit- 


Wigner form of the vector meson contributions in (31) with an energy-dependent 
total width for the broad p-resonance so that it vanishes at g^ < 4m^ and adopting 
constant total widths Ty* for the remaining narrow resonances. 

To complete the ansatz for the hadronic dispersion relations, we have to specify 
the integrals over the hadronic spectral densities of excited and continuum states 
p(u,c)(^g'j Eq.(31). In the region below the open charm threshold, q^ = s < 4m|), 

apart from the two narrow charmonium resonances J/p and p{2S), only the interme¬ 
diate states with light quark-antiquark flavour content and spin-parity 1“ contribute. 
We make extensive use of the standard quark-hadron duality ansatz employed in the 
QGD sum rules IZ7I for the vector-meson channels. The integral over the hadronic 
spectral density p|f^(s) including and continuum states with the p and u: 
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quantum numbers is replaced by the spectral density calculated from OPE: 


p^^\s)9{s - Sh) ^ ^(imH 


ip) I 

fact,LO{u,d} 


s) + 


(s) j e{s - So) 


+ ^I'^'^frct,LO{s}(^)^(^ - ^~o) > {p = u,c-, s < 4ml )) (35) 

where only the LO contributions are taken into account, including the leading-order 
quark loops and weak annihilation diagrams. The indices {u, d} and {s} mean that 
only the diagrams with u,d and s quarks, respectively, are taken into account. The 
duality threshold sq — 1.5 GeV^ is chosen in accordance with the analysis of QCD 
sum rules in the light vector-meson channels. In the contribution of the intermediate 
ss hadronic states to the spectral density (the last term in (j^)) we include 

also the small 0-meson pole contribution. This is reflected by the choice of a lower 
effective threshold parameter sq = 4m^ ~ 1.0 GeV^. Taking at s = -|- ie) the 

imaginary parts of the loop function (12): 


1 / 9 9x / 

-lmg{q + ie,m ) = IH-^ 

TT \ T 


qz 


Ami 


(36) 


and of the WA contribution (18): 
1 




(37) 


we obtain 
Au) 


PlM = 0{s-so) 


-Cl + 2C2 + ^3 + “ 3 “ + 


247r2 




1 


247r3 

and a similar expression: 


+ ^G 4 + C's + y ) lmg{s,mi)fl^{s)e{s - sq) 


(38) 


G 3 + Y + ^5 + Y ) ■^Biri^) + 


(^GsT ^G4 + G5 + y^ Img(s,m2)/+^(s)0(s - So). (39) 


247r3 


The two above expressions specify the adopted ansatz (35) at s/^ < s < Ami- 
Note that in the LO approximation, the spectral densities (s) are real functions. 
Following Ref. [9] we slightly modify the denominator in the dispersion integral over 
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Sq < s < 4m|) replacing s — — ie ^ s — — iy/sVf,^{s,q^) where refr(s,g^) = 

7 -\/s 0 (g^ — 4m^), 7 = 0.2 is the effective energy-dependent width, where the 9 
function ensures that this width is absent at negative q^. This modihcation allows 
one to transform the smooth duality-driven spectral density towards more realistic 
series of equidistant vector mesons (cf. the model for the pion timelike form factor 
used in Ref. [2E])- The addition of NLO corrections to the LO approximation for 
the duality ansatz remains a difficult task for a future improvement, involving a 
calculation of the spectral densities of the diagrams in Figs. [^and[^ 

The spectral densities above the open charm threshold, s > 4m|), contain 

a complicated overlap of broad charmonium resonances and open-charm states, to¬ 
gether with the light-quark contributions. Moreover, starting from s = {niB 
the on-shell intermediate states with b flavour related to the imaginary part of the 
R —)■ TT form factor in the timelike region also contribute. Hence, a duality-based 
parameterization of the s > 4mj^ part of the integral over (s) will not adequately 
reflect the complicated resonance-continuum structure of the hadronic spectral den¬ 
sity. On the other hand, we only need this part of the integral at relatively small 
hence following Ref. [S], we use a simple expansion in the powers of 
truncating it at the first order: 


ds 


P 


{p)i 


4m|, 


{s - ql){s - q^ - ie) 


~ a, 


+ h 




p = u,c, 


(40) 


where au,c = and are two unknown complex parameters. 

Note that in Ref. [9] other parameterizations of the dispersion integral were also 


probed, and the results in the large recoil region were numerically close to the ones 
obtained with Eq. (40), hence we will only consider this choice. 

Finally, the dispersion relations (31) take the following form: 






?o)[ 

V=p,uj,J/i>,i>(2S) ^ 


- ql){m\, - q^ - imyF^*) 



so(so) 


ds 


_Pio(g)_ 

(s - ql){s - g 2 - iy^Feff(s)) 


\ap\ exp(i0a) \bp\ exp(i0b) 


,2 


4m|, 


(41) 


These two relation aX p = u and p = c are then separately htted to the OPE result 
obtained for the l.h.s. at < 0. After that we can use the dispersion form of 
n^p\q^) in g^ > 0 and calculate the correction to the Wilson coefficient AC' 9 (g^) 
dehned in ([^ in the whole large recoil region which we specify as: 

< g^ < . (42) 


The resulting plots are presented in Figs. and for —)■ and —)■ 

respectively. Instead of showing the £t results for directly, we present 
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Figure 9: The real part (upper left, blue online) and imaginary part (upper right, 
red online) of ACg'^{q‘^) for B~ —)■ . The solid line is the dispersion relation 

fitted to the results calculated atq‘^<0 at the central input. The shaded area indicates 
the estimated 68% C.L. uncertainties obtained in the fit to the “data” points (dots,red 
online) at q^ < 0. The values of ACg'^ {q^) averaged over the q^-bins are also shown. 
The lower panel contains the same plots for —)■ . 


the directly related, but physically more relevant plots for ACg^^\q^). At q^ above 
the .J/ijj region our approach ceases to work, mostly because the contribution of the 
hadronic dispersion integral (40) to the dispersion relation increases and the simple 
polynomial parametrization cannot be used. This is also reflected by the growth of 
the uncertainties. 

A few comments on the £t procedure of Eq. (41) are in order. Here we only 
discuss the B^ —)■ vr^decay, as the case of the neutral B decays is very similar. 
We represent the results of our calculation at negative values as “data” points 
and perform a fit for our “model function” (dispersion relation), which is valid in 
both positive and negative regions, to the obtained points. As a technical remark: 
the £t is performed by collecting in the “data” all parts of Eq. (41) which contains 
no fit parameters, i.e. only the resonance contributions (amplitudes and phases) 
and the polynomial continuum parameters are included in the “model function”. 
Furthermore, we include the error correlation of the respective points at negative g^, 
and, in addition, also of the parameter — 0)- Tfi® errors and the correlation 
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<7-(GeV^) q^(GeV^) 


Figure 10: The same as in fig. [^/or B° —)■ (upper panel) and for B° —)■ 

(lower panel). 


coefficients of the “data” are obtained by varying the input parameters within their 
intervals given in Table [T| while assuming no correlation between the parameters 
themselves. Thus, we include two correlated “data” sets in the fit for both charged 
5-meson transitions, namely, the real and imaginary parts of and We 

assume a gaussian error interval of the input parameters for this procedure and a 
maximum error correlation of 80% for the numerical stability, providing also a more 
conservative estimate. As expected, the error correlation between the “data” points 
is very large and usually exceeds 80%. In addition, we hnd a positive correlation of, 
respectively, ~ 1%(10%) for the real part and of ~ 40%(1%) for the imaginary part 
of 5*^“^ with = 0). The global minima are acceptable with Xmin = 1-93 

and Xmin — 2.53 for u and c, respectively. The central values quoted here belong to 
the global minimum, whereas the 68 % C.L. error estimate includes all minima in 
the Sx"^ < 1 region. 

7 Observables in B 

Having calculated the nonlocal amplitudes in a form of the function ACg^^\q‘^), 
we substitute this function in the amplitude ([^ of the B —)■ vrf'+f'" decay and predict 
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the observables in the accessible dilepton mass region (42). 


The only element in the complete decay amplitude, that was not specihed so far, 
is the ratio (|^ of tensor and vector 5 —)■ tt form factors entering the contribution of 
the operator. To obtain it, we evaluate the ratio of LCSR’s for both form factors 
obtained in Ref. [28]. The g^-dependence turns out to be negligible in the whole 


region of validity of the sum rules, which covers the region (42), and we obtain: 

rT(g^) ~ rT(0) = 0.98 ± 0.02. (43) 

The observables in R —>■ include the differential branching fraction, direct 

CP-asymmetry and isospin asymmetry. Note that in SM the angular distribution 
in P —)■ at hxed is reduced to an overall factor (1 — cos^ 0) in the double 

differential distribution where 0 is the angle between the momentum of the lepton l~ 
and the momentum of the P-meson in the dilepton center mass frame. In particular, 
the forward-backward asymmetry in P —)• vanishes in the SM. Hence, it 

is sufficient to calculate the dilepton invariant mass distribution of the branching 
fraction: 


1 dP(p- ^ 7r-£+r) GlalJXt 


tb- 


dq"^ 


15367r^m^ 




X • 


C, + ^Ct{q^) + 


2mh 


mg + m. 




+ |C^: 


101 


(44) 


For —>■ the corresonding formula contains r^o and an additional factor 

1/2 reflecting the normalization of the P° —)■ 7r° form factor. The resulting plots are 
presented in Figs. 11, 12 Averaging the above distribution over q\ < q^ < q^ yields 
the binned branching fraction dehned, e.g., for B~ —)■ as: 




B{B- ^7r-i+r[qlql])^ 


- ql 


dq 


^dB{B- 


—)■ TT' 


T+r) 




dq"^ 


(45) 


The predicted binned branching fractions within the region (42) are presented in 
Table 1^ for all four flavour/charge combinations. 

The most interesting characteristics of the P —)■ decay in SM is the 

dependent direct CP-asymmetry dehned for the charged P-meson modes as: 




dP(p- ^ TX-tl-)ldq^ - dP(P+ ^ Ti+tr)/dq^ 
dB{B- n-i+r)/dq^ + dB{B+ TT+i+i-)/dq^ 


(46) 


The asymmetry for the neutral P-meson modes denoted as A^p (g^) has the same 
expression with B~ —)■ P°, B~^ —)■ P°. The results obtained for this observable are 


presented in Fig. 13 
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g-(GeV^) ?^(GeV^) 


Figure 11: Dilepton invariant mass spectrum and binned branching fraction (in 
GeV~‘^) for B~ —)■ (left panel) and —)■ (right panel) with 68% 

C.L. errors (shaded region and error bars). 



Figure 12: The same as in fig. 11 for B^ 
(right panel). 


—)■ (left panel) and B^ —)■ 


Anticipating the future measurements of the g^-averaged bins of CP-asymmetry, 
we also calculate 




B(B- ^ n-i+i-[gl,qi]) - B(B+ ^ Tr+t+rlglqi]) 
B(B- ^ 7r-e*e-\qlrS) + B(B+ ^ ' 


( 47 ) 


and the analogous binned asymmetry neutral P-meson modes. 

Our predictions are collected in Table 

Finally, an important indicator of the spectator-dependent nonlocal effects, such 
as weak annihilation, is a nonvanishing differential isospin asymmetry dehned as: 




2dF(P° ^ Ti^tr)/dq^ - dT{B- TT-ti-)ldq^ 
2dF(P° ^ 7r°£+r)/dg2 + dF(p- ^ n-tt)/dq^ ’ 


(48) 


where the differential widths are understood as the CP-averaged ones. Our result is 
presented in Fig and the corresponding g^-bins of isospin asymmetry: 


M<il<il] 


2F(P° ^ 7r°£+r[g^,gi]) - F(p- ^ vr'^+r [g^, g|]) 
2F(P0 ^ nH+l-[ql ql]) + F(P- ^ 7r-£+£-[g2, g^]) 


(49) 
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(7^(GeV^) 



<7^(GeV^) 

Figure 13: Direct CP-asymmetry in —>■ (upper panel) and in —)■ 

(lower panel). 

are given in Table 

Concluding the analysis of observables in i? —)■ 7ri^i~ , we notice that the magni¬ 
tude of the predicted direct CP-asymmetry for the charged B decay modes is quite 
visible; for the neutral B decays this effect is expected to be small. In this and other 
observables our analysis generates large uncertainties in the region adjacent to J/t), 
whereas the uncertainties in the p and u: region are signihcantly smaller. This is 
partly caused by the use of QCDF to £x the relative phase between the nonleptonic 
amplitudes with p and cu, which probably leads to a slight underestimate of the 
errors in the resonance region. 


8 The B —> Tiuv decay 

The semileptonic FCNC decay B — ttuV is closely related to the charged lepton 
channel. Theoretically, this process is a very clean test of the SM, involving a 
single effective operator similar to Oio, whereas the nonlocal effects studied above 
are absent. Hence, we are in a position to predict the branching fraction of this 
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Bin [GeV2] 

[0.05,2.0] 

[2.0,4.0] 

[4.0,6.0] 

[6.0, 8.0] 

[1.0, 6.0] 

B{B-) 

'^•-‘-‘'3-0.018 

0 114+OOO8 
U.il^_0.007 

0 114+0-016 
U.ii4_0.007 

0 1 07+0-036 
-0.009 

0.126l°:Oi3 

B{B+) 

u.^^a_o 020 

0 1 rifi+0 009 

U.iOD_o 008 

0 14Q+0-016 

u.ioy_0 oil 

n 1 9Q+0.030 

U.iZO_Q Q23 

0-1681°:°]® 

2 X B{B^) 

0 140+0-009 
U.i^U_o 009 

nil 7+0-008 
'^‘-'--'-'-0.008 

u.iuy_o 008 

u.uyy_o 007 

0 iiq+0-008 

u.iiy-o 008 

2 X B{B^) 

U.iZ-9:_o 008 

U.iZ4:_0 008 

U.iiO_o 007 

u.iuy_o 008 

n 1 91 +0-008 
W.1Z1_0 008 

•^CP 

— 0 171+0-027 
'J--'-'-‘■-0.045 

-0.1561°;°^ 

-0.099tng 

U.uyi_o 053 

-0.143l°;°i 

Aoo) 

•^CP 

u.uoa_o 015 

-0.0281^}° 

-0.028to+5 

1 

0 

0 

1 + 

0 0 

b b 
to to 

CO CO 

-0.008+°}° 

Ai 

n 1 qc:+0.033 
U.iyo_0 035 

-0.0201°;°31 

-0-02+ni 

—0 n2T'"°'°®*^ 
U.UZi-0 100 

-0.0631°;°°° 


Table 4: Binned branching fractions (in units of 10 ® (GeV ^)), direct CP- 
asymmetry and isospin asymmetry of B ^ . 



Figure 14: Differential isospin asymmetry for B —)■ decays. 


decay with a better accuracy than for the B —>■ . The only hadronic input 

in i? —'Kvv is the vector i? —)• tt form factor. The LCSR form factor [19] given in 


(27), provides an extrapolation beyond the large recoil region up to the kinematical 


limit = {ms — revealing a good agreement with the lattice QCD results in 

the low recoil region. We use this form factor to predict the total branching fraction 
of the B —)• TTuu decay. 

The effective Hamiltonian encompassing the b —)■ dui? transition in the SM can 
be written as: 


njb^di/u 

^eff 


AG 


V2 


F \ ^ ®em 




An 




with the (scale-independent) Wilson coefficient: 

1 


CiOv — 


sin^ Qw 


Xo{xt) + —Xi{xt)^ , 


(50) 


(51) 
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Figure 15: The differential branching fraction of the B tt uu decay 


where Xt = ml/rri^ and the functions Xo(x) and Xi[x) can be found in Ref. [12]. 

The differential branching fraction of the B~ —)■ decay summed over neu¬ 

trino flavours has the form: 


1 dB{B —)■ 71 i/v) 
tb- dq^ 


1 

tb- 


\ dB(yB —^ 

£=eyfl^T 
2 


dq^ 


(52) 


2567r^m^ 




Substituting the form factor f^-K^T) from Eq. (27), the numerical values of the 
Wilson coefficient = —6.79 and other parameters in Eq. (52), we obtain the 
differential branching fraction shown in Fig. 15 Integrating it over 0 < < (m^ — 

m.,rY we obtain: 


B{B- TT-uu) = 2B{B^ 7r°z/P) = (2.39l|]:^^) x 10“^. (53) 


Despite the fact that this branching fraction is well within the reach of R-physics 
experiments, a severe problem is the identihcation of the hnal state with respect to 
the background. 


9 Conclusion 

In this paper we calculated the hadronic input for the rare FCNC decay B —)■ 
in the large recoil region of the pion, i.e., at small and intermediate lepton- 
pair masses up to the J/'0 mass. We focused on the most difficult problem in the 
theory of these decays: the effects generated by a nonlocal overlap of the pointlike 
weak transition with the electromagnetic lepton-pair emission. At > 0 the nonlo¬ 
cality involves long distances, including the formation of hadronic resonances - the 
vector mesons. On the other hand, this part of the decay amplitude is not simply a 
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background for the FCNC b —)• di~^£~ transition, but provides the strong-interaction 
phase. The latter, combined with the CKM phase, generates the unique characteris¬ 
tics of the B —)■ decay in SM, that is the g^-dependent direct CP-asymmetry, 

suppressed in the b —)■ si^i~ decays. 

To avoid the complications related to the long-distance part of the nonlocal 
effects, we employed the method used earlier in Ref. [H] for B —)■ Ki^i~. The nonlo¬ 
cal contributions to i? —)■ transitions were calculated one by one, combining 

QCDF and LCSRs at spacelike < 0, where the quark-level diagrams are well 
dehned and the nonlocality is effectively reduced to the distances of 0{1/^/\^\). 
We also used the recently updated [19] R —)■ vr form factor from LCSR. The accu¬ 
racy of our calculation is characterized by taking into account, in addition to the 
factorizable quark-loop effects and the factorizable NLO corrections, also the impor¬ 
tant nonfactorizable contributions: the soft gluon emission, spectator scattering and 
weak annihilation. We then combined the quark-level calculation with the hadronic 
dispersion relation and htted the parameters of the latter to access the > 0 region. 
The main result of our calculation is presented in a form of the g^-dependent and 
process-specihc correction to the Wilson coefficient of the semileptonic 

operator Og. Apart from the numerical prediction for /S.C'^^^\q^), we also estimated 
the uncertainties due to the input parameter variation. We predicted the observables 
in R —)■ tt£'^£~ , including the differential branching fraction, direct CR-asymmetry 
and the isospin asymmetry. The main advantage of the method used in this paper is 
the possibility to access the p, u resonance region and, simultaneously, to approach 
the charmonium region from below. 

The accuracy of the calculation carried out in this paper can be improved further. 
On the theory side it is worth to calculate the nonlocal contributions using entirely 
LCSRs instead of the QCDF approximation. This will allow one to assess the missing 
power corrections. Such analysis is possible at least for the weak annihilation and 
for the hard spectator contributions. A more elaborated ansatz for the hadronic 
dispersion relation, including the radial excitations of light vector mesons, is also 
desirable. For that, more accurate data on the R —)■ Rtt nonleptonic decays and a 
better understanding of the structure of various nonleptonic amplitudes are needed. 

Let us compare our results with the two most recent analyses of the R —>■ 7i£^£~ 
decay. In [22], only the factorizable nonlocal contributions were taken into account, 
approximated by the quark-level diagrams at positive g^, embedded in the short- 
distance coefficients. Only the differential branching fraction was calculated, with no 
prediction for the OR-asymmetry. In [6], the QCDF method was systematically used 
at positive g^, therefore the resonance region of g^ was not accessible. In the region 
between 2 GeV^ and 6 GeV^ the branching fraction obtained in [6] is somewhat 
smaller than our result, whereas the OR-asymmetry is close to our prediction. 

We emphasize that our method produces a quantitative estimate of the nonlocal 
effects in the whole large-recoil region, starting from the kinematical threshold of 
the lepton-pair production. The price to pay is a model dependence of the ansatz for 
the dispersion relation, related to the nonleptonic R —)■ Rtt decays. The function 
ACg^^\q‘^) obtained in this paper can be used in further analyses of the R —)■ 
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decay, e.g., when adding to the decay amplitude in SM certain new physics 
contributions. But hrst of all, it will be very interesting to confront our prediction for 
the direct CP-asymmetry in i? —)■ with the data. Note that the b —)■ d£^i~ 

effective interaction is also probed in —)■ decay. Its branching fraction 

measurement by LHCb and CMS collaborations [30] still leaves some room for new 
physics, making further studies of 6 —)■ decays very important. 
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Appendix A: Operators and CKM parameters 

In Table we list the operators entering the effective Hamiltonian ([^, and their 
Wilson coefficients calculated at LO for three different renormalization scales, where 
ttem = e^/(47r) is the electromagnetic coupling, gs is the strong coupling. We use 
the standard conventions for the operators Of (p = u, c) except the labeling of Of 
and Of is interchanged, as in jH). In the quark-penguin operators q = u,d, s, c, b and 
the mass of d-quark in Oj-y and Osg is neglected. The sign conventions for covariant 
derivatives, q-matrices, left- and right-handed components of the quark helds are 
the same as quoted in the Appendix of [9]. The electroweak parameters used to 
calculate the coefficients Ct are [18] 

ttem = — , sin2(0vi.) = 0.23126 , mw = 80.385 GeV, 
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Gf = 1.1663787 x lO"® GeV^ , m, = 91.186 GeV, m* = 173.3 GeV. (54) 

We use the GKM mixing matrix in term of Wolfenstein parameters 


UCKM 


^ 1 - AV2 

-A 

AA^(1 — p — irj) 


A A\^{p — ip) ^ 
1 - AV2 AA^ 

-AA^ 1 , 
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Operator 

p (GeV) 

2.5 

3.0 

4.5 

= {dLlpPl) (pLl^bL) 

Oi 

1.169 

1.148 

1.111 

O 2 = {dilpp’L) 

02 

-0.360 

-0.324 

-0.255 

O 3 = (dljM E 
<? 

C'3(x10-2) 

1.700 

1.503 

1.144 

^>4 = {dhvK) E 
<? 

C'4(x10-2) 

-3.602 

-3.271 

-2.630 

O 5 = {dLlpfL) E {PrI^Qr) 

C'5(x10-2) 

0.985 

0.910 

0.756 

Oe = {dijpbf) E 
q 

C'6(x10-2) 

-4.829 

-4.258 

-3.236 

= -10 {dLa>^''bR) Fp, 

cf 

-0.356 

-0.343 

-0.316 

0 >s = -Wt G”"" 

r^eS 

-0.166 

-0.160 

-0.150 

0 , = ^ {dr^rbr) ih/) 

Og 

4.514 

4.462 

4.293 

Oio = ^ {dL^bL) 

Olo 

-4.493 

-4.493 

-4.493 


Table 5: Effective operators and Wilson coefficients. 


taking into account that p ~ p (1 + A^/2) and p ~ p (1 + A^/2) and using the current 
values [IH] obtained from the global CKM £t: 

A = 0.22537 ± 0.00061, A = 0.814;°;°23^ 

p = 0.117 ±0.021, p = 0.353 ±0.013. (55) 

This results in the following combinations of CKM elements we use: 

A^/At = -0.0274 - i 0.3896, \\u/\t\ = 0.3906, arg(A„/At) = -94.02° 

A^/Ai = -0.9719 ± i 0.3998, |Ac/Ai| = 1.0509, arg(Ac/At) = 157.64° . (56) 

Appendix B: Amplitudes of B > p{uj)7T in QCDF 

Here we present the expressions of the QCDF amplitudes uni for the B —)■ 
(p°, u)7r~ nonleptonic decays. Our operators differs from the ones in [10] by a factor 
1/4 whereas the labeling of Oi ^2 is the same. The expressions for the parts of 
B~ —)■ p°7r“ and B~ —)■ cutt” amplitudes multiplying A^ (p = u, c) are: 

^ b - p - k - = Al^p(^^P«[«2(vrp) - /32(vrp)] - a^(7rp) - /3f(7rp)^ 

+Ap^ ( Spu [ai(pvr) ± /32(p7r)] ± ^^(pvr) ± (pvr) j, (57) 
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^5pn[a2(vra;) + /?2(7ra;)] + 2a^(7ru;) + a^(7ruj) + /3f(7ra;)^ 

+A^^^5p„[ai(a;7r) + /32(cu7r)] + a^(uj7r) + /?f(a;7r)^ , (58) 

where the factorized combinations of form factors and decay constants are: 

ATrp(uj) = B'Ki^p)fpiA-: ^p{io)TT = fn ■ ( 59 ) 

In addition to the already introduced notation, ^sp(cj) (0) fh® above is the relevant 
B —)■ p{u) form factor taken at = 0, neglecting the pion mass squared; for the 
B ^ TT form factor we approximate The parameters af(MiM 2 ) are 

defined as follows [TO] : 


ai{MiM2) — ai{MiM2), i — 1,2, 




a. = 


a^(MiM 2 ) + af (M 1 M 2 ), if M 2 = p, cj, 
al{MiM 2 ) - af (M 1 M 2 ), if M 2 = vr, 

<(MiM2) + r^^al{M,M2), if M 2 = p, (u, 
<(MiM 2 ) - rf 2 aP(MiM 2 ), if M 2 = vr. 


where 


2ml 


= 






p,ljJ J p,iJ 


^ mb{mu + md)' ^ mb fp,uj 

and is the vector-meson transverse decay constant, defined as 

{0\q^^’'q\V{q))=tk^{eW-^v<infv 

with fc-*- = 1/a/2 ioT q = u, V = p^,u. 

The quantities af(MiM 2 ) have the form [2|: 


a\{MiM2) — [Ci + 


iVc 

Ci±i CpOis 


iV,(M2) 


+ pnMi) . 


(60) 

(61) 

(62) 

(63) 

(64) 

(65) 


where the upper (lower) signs apply when i is odd (even); iVj(M 2 ) = 0 for i = 6 and 
M 2 = p, to and Ni{M 2 ) = 1 in all other cases. The parameters /3f(MiM2) involve 
the weak annihilation contributions: 


/?f(MiM2) 


-/b/mi/m2 


(M 1 M 2 ). 


( 66 ) 
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The expressions used for the separate contributions in Eqs. (65), (66): Vi{M 2 ) 
(one-loop vertex correction), Hi{MiM 2 ) (hard-spectator scattering), Pf{MiM 2 ) (pen¬ 
guin contractions) and {M 1 M 2 ) (weak annihilation) can be found in OITO]. They 
were calculated in QCDF in terms of the perturbative kernels convoluted with the 
DAs of the B meson, pion and p{u:) meson. The latter DAs include the Gegenbauer 
moments and similar to the ones that are contained in the pion twist-2 
DA (|^. 

For the numerical analysis of the B~ —)■ amplitudes we need additional 

input parameters listed in Table where, in order to decrease the uncertainty, the 
^^“^(0) form factor is calculated multiplying the ratio A^^(0)//|^’^^(0) obtained 
from the LCSRs with the i?-meson DAs [23] with the form factor taken 

from the most accurate LCSR with pion DAs |19j . 


Parameter 

Ref. 

{mu + mrf)(lGeV) = 7.0^^;^ MeV 

Iig 

= 0.316 ±0.021 

ng 

= 0.396«;“|; 

<“(0) ~ ^='■(0) 

ra 

f^ilGeV) = (0.160 ±0.010) GeV 
/j((lGeV) = (0.145 ±0.010) GeV 
a^’‘^(lGeV) = a^’"^^(lGeV) = 0.09l|];J? 

1211 


Table 6: Additional input parameters related to the light vector mesons and used in 
the QCDF amplitudes of B ^ decays. 
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